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1. Introduction 

Idealisations play an important role in modelling a wide range of physical situations. 
In many field theories a particularly useful idealisation is to replace an extended 
source which is concentrated in a small region of space by a point charge. Such an 
approximation is physically reasonable provided the internal structure of the source can 
be neglected. In a similar way sources concentrated near a line or a surface can be 
described in terms of strings or shells of matter. On trying to describe this idealisation 
mathematically it is natural to use a delta function to describe the source, and hence in 
this context distributions arise in a natural way. Distributions are also used to describe 
a number of other important physical scenarios such as the description of shock waves 
and the junction conditions between matter and vacuum regions. 

In the case of linear field theories such as electrodynamics, distribution theory in 
fact furnishes a consistent framework which has the following two important features. 
Firstly, since the Maxwell equations are linear with respect to both sources and fields, 
one can allow distributional solutions of the field equations as well as classical smooth 
solutions. Secondly it is guaranteed that smooth charge densities that are close to those 
of a point charge (in the sense of distributional convergence) produce fields that are close 
to the Coulomb field. While the first property permits one to have a mathematically 
sound formulation of concentrated sources it is precisely the latter notion of "limit 
consistency" which renders the idealisation physically reasonable. 

One would like to have a similar mathematical description of concentrated sources 
in the theory of general relativity. However, general relativity is different from other field 
theories in two important respects. Firstly the Einstein field equations are nonlinear, 
so that one cannot simply pass from smooth solutions of the field equations to weak 
solutions as one can in a linear theory. More precisely the curvature tensor is a nonlinear 
function of the metric and its first two derivatives. Thus the metric needs to be at least 
to guarantee that the curvature is continuous. Mathematically one can weaken this 
condition to allow a (i.e. first derivative Lipschitz continuous) metric and this is 
sufficient for most results in differential geometry to remain valid. In certain special 
situations one can lower the differentiability requirements still further and formulate 
the field equations in a way which avoids ill-defined products of distributions. For 
example it is possible to describe shells of matter [32], JH] and gravitational radiation 
inn] within the context of classical distribution theory. Using the standard definition of 
the curvature 

pa -pa -pa , pa pf pa pf ( -i \ 

^ bed — ^ db,c ^ cb,d cf^ db ^ df^ cb \^) 

one sees that one wants the connection to have a (weak) derivative and be locally square 
integrable in order for the left hand side to make sense as a distribution. This led Geroch 
and Traschen j^H] to introduce a class of "regular metrics" for which the components 
of the curvature are well defined as a distribution. Geroch and Traschen went on to 
show that such regular metrics can only have curvature with singular support on a 
submanifold of co-dimension at most one. Thus for a 4-dimensional spacetime a metric 
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representing a shell of matter could belong to this class but a string or particle could not. 
We will consider such distributional solutions of Einstein's equations in more detail in 
§2. The second important respect in which general relativity is different from other field 
theories is that one does not have a fixed background metric, but instead the geometry 
is determined by the field equations. Spacetime is described by a manifold together 
with a Lorentz metric which is assumed to be sufficiently different iable for Einstein's 
equations to be defined. One then detects the presence of singularities by showing that 
the spacetime is incomplete in some sense. The problem with this approach is that many 
physically reasonable spacetimes contain singularities according to this definition. For 
example solutions to Einstein's equations representing cosmic strings are singular. What 
one would like to do is to lower the differentiability required of the metric to permit a 
wide class of "distributional geometries" which represent physically reasonable solutions 
and are also mathematically tractable. Owing to the nonlinear nature of Einstein's 
equations such distributional geometries will in general require some nonlinear theory 
of generahsed functions. Moreover, as pointed out by Isham [22] distributional solutions 
are not only important in allowing one to deal with concentrated sources or describing 
weak singularities but (based on our experience with linear field theories) should also 
be included in any path integral formulation of quantum gravity. 

Although there has long been a desire to allow distributional geometries in general 
relativity it is only comparatively recently that any real progress has been made in 
realising this aim. Many of the early attempts (e.g. IHHl, [IS], 0) used methods 
which were specifically adapted to the particular problem being considered and whose 
general applicability was uncertain. More recently a number of different authors (for 
an overview see [22], jlHl Ch. 5.2]) have investigated distributional geometries using an 
approach based on Colombeau algebras. These were developed by J.F. Colombeau in the 
1980 's and contain the smooth functions as a subalgebra and the distributions as a linear 
subspace. The key idea of this approach is regularisation through smoothing and the 
use of asymptotic estimates with respect to a regularisation parameter. In particular, 
it provides a mathematically consistent way of multiplying distributions and a unified 
view on calculations involving various regularisation procedures. An important feature 
of these algebras is the notion of association which gives a correspondence between 
elements of the algebra and distributions. This allows one to use the power of the 
algebras to do mathematical calculations but then use the concept of association to 
interpret the final result in terms of classical distributions and give the solutions a 
physical interpretation. 

In §2 we review the extent to which it is possible to incorporate nonlinear operations 
into classical distribution theory. We will show that some very limited operations 
are permitted and that some apparently reasonable "multiplication rules" lead to 
inconsistencies. We then consider the implications of this for distributional solutions of 
Einstein's equations and end the section by looking in detail at the properties of the 
Geroch-Traschen regular metrics. In §3 we give a brief review of Colombeau theory and 
explain how it is able to circumvent the result of Laurent Schwartz on the impossibility 
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of multiplying distributions. Historically one of the first singular solutions of Einstein's 
equations to be studied from a distributional viewpoint was the Schwarzschild solution 
which was considered by Parker in 1979 [HH]- The key observation is that when written 
in Kerr-Schild coordinates the coefficients of the metric are locally integrable. The 
Schwarzschild and Kerr solutions have been studied in these coordinates by Balasin 
and Nachbagauer in a number of papers and in §4 we describe their work using the 
language of the (special) Colombeau algebra. If one boosts the Schwarzschild solution 
with velocity v then one obtains, in the limit that the velocity is that of light, the 
"ultrarelativistic Schwarzschild solution". This metric was investigated by Aichelburg 
and Sexl pp who showed that it could be thought of as a gravitational shock wave. When 
one does a similar calculation with the ultrarelativistic limit of the Reissner-Nordstr0m 
solution one obtains a solution with vanishing electromagnetic field but (5-function energy 
density. In §5 we show how Steinbauer was able to explain this physically surprising 
result using Colombeau algebras. In §6 we consider geodesic equations for impulsive 
gravitational wave spacetimes and show how Colombeau algebras provide an appropriate 
way of obtaining the results expected on physical ground without the need to make use of 
ad hoc "rules" for the multiplication of distributions. The study of conical singularities 
is another area where various authors had used a variety of regularisation procedures to 
obtain the physically plausible result that a cone has 5-function curvature at the vertex. 
However this result appeared to be at odds with the results of Geroch and Traschen. In 
§7 we review the analysis and resolution of the problem by Clarke, Vickers and Wilson 
using the (full) Colombeau algebra. In §8 we discuss questions of coordinate invariance 
of Colombeau algebras and review work on global and diffeomorphism invariant versions 
of the construction leading to a "nonlinear distributional geometry" . In §9 we look at 
"generalised hyperbolicity" and how weak singularities may be regarded as distributional 
solutions of Einstein's equations. Finally we give some conclusions and an outlook to 
future work in §10. 

2. Classical Distribution Theory and General Relativity 

In this section we briefiy review the fundamental problems encountered when one tries 
to incorporate nonlinear operations into classical distribution theory. We will then 
examine the inherent limitations this imposes on distributional products and review the 
consequences for distributional solutions of Einstein's equations. 

There has been a long history of using generalised function ideas in physics to 
model point sources and discontinuous phenomena. Such generalised functions were 
put on a sound mathematical footing by the development of the theory of distributions 
through the work of S. L. Sobolev [ZH] and L. Schwartz jUj. The basic idea is to make 
distributions dual to a space of smooth "test functions". To introduce some notation 
we let T>(R"') denote the space of smooth functions of compact support on M". If 5* is 
a linear form S : ©(R") —>■ C then we will denote the action of S* on G P by {S, (p). 
The vector space of distributions V'iW") is then defined to be the space of continuous 
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linear forms on X>(R"). In a similar way distributions on an orientable manifold M are 
defined as continuous linear functionals on the space of compactly supported n-forms 
i.e. V{M) := [r2"(M)]'. A rich theory of distributional tensor fields (and sections of 
more general vector bundles) has been developed by Marsden in [^21; for a pedagogical 
introduction see [2H1 Ch. 3.1]. 

The theory of distributions quickly proved to be extremely successful both in 
applications to the study of linear partial differential equations and in justifying the 
use of generalised functions in physics. For example the Malgrange-Ehrenpreis theorem 
shows that any linear PDE with constant coefficients has a fundamental solution in the 
space of distributions. However the theory soon displayed its natural limitations when 
in 1957 H. Lewy gave his famous example of a linear partial differential equation 
with smooth (non-constant) coefficients which does not have a distributional solution. 
A second difficulty with the theory is that the definition of a distribution as a linear 
functional does not make it easy to define the product of distributions. This prevents one 
from using distributions to investigate nonlinear PDEs with singular data or coefficients. 

Although the Lewy example of a linear PDE without a distributional solution 
came as a great surprise, the difficulties that one encounters with the multiplication 
of distributions are much easier to understand and can be seen by looking at some 
simple cases. In the following we shall briefly discuss such examples concentrating on 
the powers of the Heaviside function H and its product with Dirac's delta function 
6. If we regard if as a discontinuous function then if™ = H [m E N). However if 
we differentiate this formula and use the Leibniz rule for the derivative the following 
one-line calculation implies the vanishing of the delta function: 

2H6 = {H^y = {H^y = 3H^6 = 3H6, so H6 = hence 6 = 0. (2) 

Another popular "multiplication rule" is 

H6 = ^6. (3) 

We demonstrate the problematic nature of this rule (and in fact any rule of the form 
HS = c6, c E M) by considering the simple ODE 

y\t) = 6{t)y{t), y(-oo) = 1, (4) 

(for an amusing discussion of this equation see |SI1)- 

Using the ansatz y(t) = 1 + aH{t) and Q we flnd a5 = {1 + a/2)5, hence a = 2 
and the solution takes the form 

y{t) = l + 2H{t). (5) 

On the other hand a different approach motivated by the requirement of stability 
under perturbation is to regularise the singular coefficient by a sequence 5n weakly 
converging to 5. Then the solution to the regularised equation is given by ?/„(t) = 
exp(J^ 5n{s)ds) + 1 which converges to 

y(t) = l + {e-l)H{t), (6) 
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which obviously does not coincide with the previous solution. 

The deeper reason behind these and all other inconsistencies in this realm is the 
incompatibility of the laws of a (commutative) differential algebra with the formulae 
H"^ = H and H' ^ 0. This insight was put into stringent form by L. Schwartz himself 
in his incompatibility result j26j, which says that if the vector space V of distributions is 
embedded into a differential algebra {A, +, o) then the following properties are mutually 
contradictory: 

(i) V is linearly embedded into A and f{x) = 1 is the unity in A. 

(ii) There exist linear derivation operators di : A ^ A satisfying the Leibniz rule. 

(iii) 9j|x)' is the usual partial derivative. 

(iv) o|^fex(7fe (for k finite) is the usual pointwise product of functions. 

It was this result that led to the idea that it was impossible to multiply distributions. 
However given that repeatedly differentiating a function eventually produces a 
distribution it is perhaps unreasonable to insist on (iv) but instead we should only 
require that the product of smooth functions is the usual pointwise product. We will 
see in §3 that this is precisely the condition satisfied by Colombeau algebras. 

In the rest of this section we examine the extent to which one can apply linear 
distribution theory in the context of general relativity. After summarising work done 
in a number of special cases we review a classical paper by Geroch and Traschen [20] 
in which they set up a "maximal" distributional framework by finding the "largest 
possible" class of spacetime metrics which allow for a distributional formulation of the 
field equations, and we discuss its limitations. 

Spacetimes involving an energy-momentum tensor supported on a hypersurface of 
spacetime have long been used in general relativity (see [SDl EH 120] and IHEl lEE], as 
well as the references therein). Consider a submanifold 5* of codimension one dividing 
spacetime into a "lower" and "upper" part and let the metric be smooth up to and 
including S from each of its sides but allow for a jump of its first derivatives across 5*. 
Writing out the Einstein equations in terms of the extrinsic curvature of S one finds 
junction conditions very similar to the ones in electrodynamics (see e.g. [Z^, Ch. 3.7). 
More precisely, the jump of the extrinsic curvature is interpreted as the surface stress- 
energy of a surface layer located at S. In the case of S being timelike this arrangement 
represents a thin shell of matter, while if S is null it may be interpreted as a thin shell 
of radiation (see e.g. [Snj). In jSH] W. Israel has given a general formulation of this 
widely applied approach with the practical advantage that no reference to any special 
coordinate system is required; the four- dimensional coordinates may be chosen freely 
and hence may be adapted to possibly different symmetries in the upper and lower part 
of spacetime. 

On the other hand Lichnerowicz jSS], has given an alternative description using 
tensor distributions assuming the existence of an admissible continuous coordinate 
system across S. This formalism was used by Lichnerowicz |S3J EH] and Choquet- 
Bruhat to study gravitational shock waves. They derived algebraic conditions on 
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the metric across the shock (the "gravitational Rankine-Hugoniot conditions") as well 
as equations governing the propagation of the discontinuities. The respective formalisms 
of Israel and Lichnerowicz were shown to be equivalent in 

The description of gravitational sources supported on two-dimensional submanifolds 
of spacetime, however, is more delicate. Israel [HZj has given conditions under which a 
sensible treatment of the field of a "thin massive wire" is possible. He isolated a class of 
"simple line sources" which possess a linear energy-momentum tensor and hence allow 
a well-defined limit as the wire's radius shrinks to zero. We will look at line sources in 
more detail below. 

On the other hand Taub |H1] has claimed to have generalised Lichnerowicz's 
formalism to include gravitational sources supported on submanifolds of arbitrary 
codimension in spacetime. However, he had to fix the ill-defined products by 
"multiplication rules", in particular, by using equation (jHl). 

We now begin to review the systematic approach by Geroch and Traschen in 
analysing the structure of the nonlinearities of the field equations to see how far one 
can get avoiding ill-defined distributional products. More precisely, the quest is for a 
class of metrics allowing for a distributional formulation of the Einstein tensor in order 
to assign to the spacetime — via the field equations — a distributional energy-momentum 
tensor representing the "concentrated" source. Note that there are two contradictory 
demands on this class of metrics: on the one hand these metrics should be "nice enough" 
to permit the distributional calculation of the curvature entities, while on the other 
hand they should be "bad enough" to have the Einstein tensor and hence the energy- 
momentum tensor concentrated on a submanifold of a high codimension in spacetime. 

We write out the coordinate formula of the Riemann curvature tensor in terms of 
the Levi-Civita connection and the connection in terms of the metric 

Rabc"^ = 2rg[^r^]^ 29[6rf]^ , (7) 

= g^'idime - ld,9bc) . (8) 

and try to "save" these equations by putting just as much restrictions on the metric 
tensor as needed to allow for a distributional interpretation of the respective right hand 
sides. For the first term in ((7j) it is obviously sufficient to assume F^^ to be locally square 
integrable. Since Lj^^ C Lj^^ this requirement actually also suffices to interpret the 
second term in ((Tj) as the weak derivative of the regular distribution F^^. Furthermore 
from equation (jH)) we see that it is sufficient to demand g"-^ to be bounded locally 
almost everywhere in order to produce locally square integrable F^^ from locally square 
integrable first weak derivatives of Qab- This motivates the following definition. 

2.1 Definition. A symmetric tensor field Qab on a four- dimensional manifold M is 
called a gt-regular metric if gab and g""^ G Lj^^ fl H^^^ . 

In the above definition Lj^^ denotes the space of locally bounded functions and H^^^ 
denotes the Sobolev space of functions which are locally square integrable and also have 
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locally square integrable first (weak) derivative. Note that although the above definition 
appears to be stronger than that in it is actually equivalent to the original one where 
it was merely demanded that Qab as well as g""^ are locally bounded almost everywhere, 
and the first weak derivatives of gab are locally square integrable. In fact, these assertions 
imply dcg"''' G L^^^ since g"-^ = cof (gab) / det gab, where cof (gab) denotes the cofactor of 
gab- We also note that the above conditions will always hold for a C^^ metric, for such 
a metric will always admit a locally bounded weak derivative. For a detailed discussion 
of the relationship between gt-regularity and some other regularity conditions in the 
context of axial and cylindrical symmetry we refer to 0|, chaps. 2.3-2.5. 

A further important remark on the notion of gt-regular metrics is in order. While 
the definition is coordinate invariant with the manifold fixed beforehand, in the case 
we are most interested in i.e. when we are dealing with a singular spacetime in general 
relativity, the situation is different. We are not given in advance a coordinate system 
that includes the singularity. So the question of whether a singular metric is gt-regular 
or not depends crucially on the choice of the differentiable structure which is imposed 
on the manifold to include the singular region. 

To see that gt-regular metrics actually allow for the distributional formulation of 
Einstein's equations we show that one can build a distributional Einstein tensor. To 
do this we need to show that the tensor product of the contravariant metric with the 
Riemann tensor makes sense as a distribution. By writing the terms involving the second 
derivative as a total derivative we may write this in the form 

g^fRabc ' = 2/^r:l[,r^, + 29[,(/^rf] J - 2(9[,/^)rf], . (9) 

Now the first term involves a product Lj^^ x Lj^^ hence stays locally integrable. The 
second term involves a weak derivative of an Lj^^-tensor field hence may be interpreted 
as a distribution and the third term is a product of two locally square integrable fields 
so is also locally integrable. 

Before discussing convergence for gt-regular metrics we briefiy introduce tensor 
distributions. Distributional sections of vector bundles and, in particular, distributional 
tensor fields can be defined as continuous linear forms on suitable test section spaces 
but are most easily viewed just as sections with distributional coefficients, that is 

v'liM) = v\M) ® r;{M), (10) 

where V'l and denote the spaces of distributional and smooth (r, s)-tensor fields 
respectively. 

We can now discuss an appropriate notion of convergence for gt-regular metrics. 
As already indicated above we would like the Einstein tensor of a sequence of metrics 
approximating a gt-regular one to approximate the Einstein tensor of the gt-regular 
metric. The natural notion of weak convergence for a sequence of locally square 
integrable tensor fields {{fJ']\','!jJn)n is convergence locally in square integral, i.e. 

(f^^'lX - (n ^ 00) iff y"(/i;;:::5:)„(/i£::f% t^:itX ^ ^ (ii) 
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for all smooth compactly supported (2s, 2r)-tensor densities 




Using this 



notion of convergence one may prove the following theorem. 
2.2 Theorem. (Convergence of gt-regular metrics) 

Let Qah and {{gab)n)n be a gt-regular metric and a sequence of gt-regular metrics 
respectively and let 

(i) {{gab)n)n o^^d {{g"'^)n)n be locally uniformly bounded, and 

(ii) {9ah)n 9ah, {9''^)n ^ 9°^, and {dagbc)n dagbc locally in square integral. 
Then {RabcX ^ Rab/ in V'\{M) and hence {Gab)n ^ Gab m V'l{M). 

We note that the space of gt-regular metrics is complete with respect to the notion 
of convergence defined by hypotheses (i) and (ii). Moreover, let gab be a continuous 
gt-regular metric then there exists a sequence of smooth metrics {{gab)n)n converging to 
gab in the sense of (i) and (ii). 

Before actually checking which class of gravitational sources may be described by 
gt-regular metrics we start with the following heuristic consideration of the behaviour 
of gt-regular metrics. Suppose S is d dimensional submanifold of a 4-dimensional 
spacetime M and the metric gab is smooth on M \ 5" but some of its components diverge 
as one approaches S. What order of divergence is allowed if gab is to be gt-regular? 
Let r be a typical distance from S measured by some background Riemannian metric 
hab and suppose the components of gab diverge at the rate of for some positive 
number s. Then the weak derivatives of gab diverge like r~^~'^ while the volume element 
is proportional to r^"*^. In order for the derivatives of the metric to be locally square 
integrable we therefore require that 2(— s — l) + Z — d> —1, and hence that 



Hence we see that the components of gt-regular metrics must grow more slowly than 
a rate of r~^'^'^^'^ as one approaches a d-dimensional submanifold in 4-dimensional 
spacetime. In particular, the larger the codimension of the submanifold the more 
strongly the components of the metric may diverge. However, as shown by the following 
theorem, there are severe constraints on the dimension of S. 

2.3 Theorem. (Concentrated sources from gt-regular metrics) 

Let S he a submanifold of dimension d = (0, 1, 2, 3) of a four- dimensional manifold M 
and let ^j^.'.'.j^ ^ a tensor distribution satisfying 



(ii) Tj^',''jl is the sum of a locally integrable tensor field and the weak derivative of a 
locally square integrable tensor field (hence is of the form of the Riemann tensor of 
a gt-regular metric). 




(12) 



(i) supp(T;^^:;j;) C S, and 



Then d = 3. 
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This theorem fits the picture described earher in this section, i.e. that gravitating 
sources with their support concentrated on a 3- dimensional submanifold have been 
treated successfully in the literature while sources concentrated on submanifolds of 
higher codimension have turned out to be more subtle to deal with. However, it 
should also be emphasised that although the gt-regularity conditions are sufficient for 
the curvature to make sense as a distribution they are certainly not necessary as shown 
by the examples below, so that the above result should not be interpreted as implying 
that a spacetime with distributional curvature can only have the source confined to a 3- 
dimensional submanifold. Indeed in certain algebraically special situations — and using 
a preferred coordinate system — some of the curvature quantities may be defined for non 
gt-regular metrics. For example impulsive pp-waves ^H] have been treated extensively 
using distributions. 

Summing up, gt-regularity provides us with a large class of "badly behaved" metrics 
which nevertheless allows one to formulate the field equations in a "stable" way. We are, 
however, sailing close to the wind as may be seen from the fact that energy conservation 
may not be formulated in general for gt-regular metrics. Indeed, the left hand side of 
the Bianchi identities V[aRbc]de involves a product of the distributional coefficients of the 
Riemann tensor with the non-smooth Christoffel symbols and this is only well defined 
if one imposes additional conditions on the metric. Similarly, gt-regular metrics may 
not be used to raise or lower the indices of a general tensor distribution since the tensor 
product would again involve a multiplication of distributions. 

By looking in more detail at the combination of terms that one has in the expression 
for the curvature, Garfinkle in |2S1 has generalised the formalism of Geroch and Traschen 
to include a slightly more general class of metrics. However, in extending the class 
of metrics in this way one can no longer establish the convergence theorem which 
one has for gt-metrics. One is therefore forced to give up the requirement of "limit 
consistency". This is another indication that with the Geroch Traschen definition of 
regularity one has gone about as far as possible using conventional distribution theory. 
Staying strictly within the mathematically and physically consistent setting given by 
this theory, one has to restrict oneself to a class of metrics that excludes physically 
interesting cases such as strings and point particles. If one wants to describe more general 
gravitational sources the nonlinearity of the field equations forces one to go beyond the 
limits of classical distribution theory and face true conceptional problems. A consistent 
framework allowing for nonlinear operations on singular (e.g. distributional) objects is 
provided by Colombeau's algebras of generalised functions, which we introduce in the 
next section. 

3. A brief review of Colombeau theory 

In this section we give a brief introduction to Colombeau algebras. For more details see 
[T7], [IH], [12], |2H|- As we said in the previous section the definition of distributions as 
linear functionals is not well suited to formulate a definition of multiplication. However 
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it is common to visualise the Dirac delta function as the limit of a sequence of smooth 
functions, all with integral one, whose support gets concentrated at the origin. In fact 
it is possible to give these ideas a precise mathematical formulation and an alternative 
sequential approach to distribution theory was developed by Mikusihski [Hlj as early as 
1948 (see also Temple ^3)- ■'■^ ^^^^ approach a distribution is an equivalence class of 
weakly converging sequences of smooth functions modulo weak zero-sequences. Working 
with a more subtle quotient construction Colombeau was able to construct a differential 
algebra A satisfying conditions (i)-(iii) of §2 but with (iv) replaced by 

(iv') o|(^ooxc°° is the usual pointwise product of smooth functions. 

In order to introduce the basic concepts we will start by describing the special (or 
simplified) Colombeau algebra on M". The basic idea is to consider generalised functions 
as 1-parameter families of smooth functions {fe}- Our basic space will thus be 



We now want to consider how we can represent a function / of finite differentiability as 
an element of this space. If we start with some $ G ^^(IR") with integral one then we 
can rescale this to obtain a family of functions 



with the property that 6 in V as e 0. Hence if we take the convolution of / 

with we obtain a family 



of smooth functions that converge to / in V as e tends to zero. We will refer to $ in 
the above expression as a molifier. 

Of course we can also apply the above formula to a smooth function /. However 
for a smooth function we can also represent / as an element of S by considering the 
constant family fe{x) = f{x). For the case of a smooth function we would like both these 
possible representations to be equivalent. Using a Taylor series expansion to compare 
the difference between these expressions we are lead to define two representations to 
be equivalent if they differ by a "negligible function" which is defined as a 1-parameter 
family of functions which on any compact set vanishes faster than any given positive 
power of e. Since we are trying to construct a differential algebra we also require the 
derivatives of / to satisfy this property and the resulting set Af to be an ideal. Clearly Af 
is not an ideal in £^(R"), but by restricting this space to "moderate functions" £^Af(IR") 
which grow no faster than some inverse power of e, one does have an ideal and we may 
define the differential algebra Q as the quotient. 




(13) 




(14) 




(15) 



3.1 Definition. 
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(i) (Moderate functions) 



£^(M") := {{/J : WK CC M", Va G N[J, 3p G N stic/i t/iat 
sup iD^'f^ix)] ^ 0{e-P) ase^O}. 



(a) ( Negligible functions ) 



Af{W) := {{/J : \/K CC M",Va G NJJ,Vg G N such that 
sup|Z^"/,(x)| ^0(e«) ase^O}. 



(Hi) (Special algebra) 



Note K CC M" indicates that K is compact and we have also employed the standard 
multi-index notation for -D"/. 

Thus a nonlinear generalised function / denoted by / = [{fe}] is an equivalence class 
of moderate sequences of smooth functions modulo negligible ones; it is represented by a 
moderate sequence of smooth functions {ft}- The space ^m(IR") is a differential algebra 
with pointwise operations and since the space of negligible functions is a differential 
ideal, Q is also a commutative differential algebra. 

The vector space of distributions is now embedded into the algebra Q by convolution 
with a molifier: More precisely we choose a molifier $ which for technical reasons 
(not to be discussed here) is a Schwartz function and has all moments vanishing i.e. 
/ $(x)x" dx = \/\a\ > 1. Then we embed T E V with compact support as 



Distributions which are not compactly supported are embedded via a localised version 
of (fTBj) using a standard sheaf theoretic construction. 

As remarked earlier one of the advantages of the Colombeau approach is that 
one may frequently interpret the results in terms of distributions using the concept 
of association or weak equivalence. A generalised function / is said to be associated to 
a distribution T G P' if for one (hence any) representative {f^} we have 



and we then write f ^ T. Note that not all elements of Q are associated to distributions. 
More generally we say for two generalised functions / ~ (7 if 



for one (hence any) pair of representatives {fe}, {Oe}- Association is an equivalence 
relation which respects addition and differentiation. It also respects multiplication by 
smooth functions but by the Schwartz impossibility results cannot respect multiplication 
in general. 



.(T) = [{T*$J]. 



(16) 




(17) 




(18) 
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The algebra presented above is the simplest of the Colombeau algebras and can 
be readily generalised to arbitrary manifolds (see §8). However it does suffer from 
the disadvantage that the embedding l of distributions and of functions of finite 
differentiability is not canonical but depends on the choice of molifier $ (see above). 
Thus one has to appeal to mathematical or physical arguments outside the theory to 
justify a particular representation of a non-smooth function. We discuss these matters 
in §8 where we also present the construction of so-called full Colombeau algebras which 
do posses a canonical embedding of distributions. In the following sections however we 
will discuss applications of Colombeau algebras to general relativity using the language 
of the special version. 



4. The Schwarzschild and Kerr Spacetimes 

In this section we review (linear and nonlinear) distributional treatments of the 
Schwarzschild and Kerr spacetimes. We use the language of the special Colombeau 
algebra although strictly speaking we should be using the special version of the theory 
of generalised tensor fields on manifolds, which we will introduce in §8. However the 
precise details will not be needed as we aim at presenting the main ideas and concepts 
in the most elementary way. 

Balasin and Nachbagauer considered rotating, charged, Kerr-Newman black-hole 
solutions in a number of papers (IH El El 13 El EH!)- The solutions considered have 
the feature that they are all examples of Kerr-Schild geometries. Such metrics may be 
written in the form 

gab = Vab + fkah (19) 

where rjab is the fiat Minkowski metric, / is an arbitrary function, ka = fjatk^ and A;" is 
null and geodetic with respect to ijab- 

The simplest example of such a solution is the Schwarzschild solution which in the 
standard Minkowski coordinates has Kerr-Schild form given by 

/ = — , F=(l,x^). (20) 
r 

The expression for the Ricci tensor of a Kerr-Schild metric takes the surprisingly simple 
form 

Rt = hV^'ri'^idedcUkdh) + dMfkdK) + dMfkeh)]. (21) 
The energy momentum tensor is then given by Einstein's equations as 

= Rl- \^IR- (22) 

One may then calculate the energy momentum tensor as follows. The Kerr-Schild form 
is regarded as being valid on the whole of Minkowski space with ^ replaced by some 
suitable regularised function. One now considers this as an element of Q and computes 
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the components of T^" (in Minkowski coordinates) as elements of Q. Finally one can 
show that 

T°^-m(5(^) (23) 

Tfe^^O otherwise. (24) 



Different authors ([HHl HOI EZl) using various regularization procedures have also 
assigned a distributional energy-momentum tensor to the Schwarzschild geometry. 
These approaches have been compared using the language of the special algebra in 

The calculation of the energy- momentum tensor for the Kerr solution is significantly 
harder. Unlike the Schwarzschild case we cannot easily write the Kerr solution as the 
limit of a one parameter family of regular Kerr-Schild metrics. The problem arises 
because of the topology of the maximal analytic extension of the Kerr solution which 
leads to a branch singularity when the metric is written using the standard "fiat" 
Kerr-Schild decomposition. Balasin [H] avoided this problem by considering metrics 
of generalised Kerr-Schild form. These are metrics which can be written 

9ab = 9ab + fkah (25) 

where gab is now a background metric, ka = Qabk^ and A;" is null and geodetic with 
respect to Qab (and also Qah because of the form of the metric). One now has 

Rl = Rl + r'r^R^cebfkdkf + \g'T[dedaUkdkb) + dbdcifk.ke) + dMfkeh)] (26) 

where -R'^^cd the curvature of cjab- 

We may write the Kerr metric in generalised Kerr-Schild form by taking the 
background metric to have the form 

Qabdx^dx^ = -d^ + dr^ + Y.de'^ + (r^ + a^)d(l)'^ (27) 

+ a"' 

where S = + cos^ 9 and (t, r, 9, (p) G x S*^. Performing all the calculations in the 
special algebra one may compute \/gRl which is found to have the following associated 
distribution 

^Rl ^ 2TT6{cos9){-a6{u)d^dub + d^d9b + m6'{u){d^ - {l/a)dl){dtb + ad(f)h). (28) 

The above calculation was carried out by Balasin using the special algebra. This 
is probably the only practical way of doing the calculation given the topology of the 
manifold and the complexity of the metric. It suffers from the usual problem when 
using the special algebra of a non-canonical embedding. Rather than embed using a 
convolution, the embedding has been chosen to preserve the generalised Kerr-Schild 
form. However the embedding used is not unique within this class and it would be 
desirable to show that any reasonable embedding which preserved the form of the 
decomposition gave the same result. It is even less clear that a "natural regularisation" in 
some other coordinate system would give the same result. Nevertheless the calculation is 
an impressive example of the complicated calculations that can be performed in general 
relativity using the special algebra. 
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5. Ultrarelativistic Black Holes 

In 1971 Aichelburg and Sexl pp derived the ultrarelativistic limit of the Schwarzschild 
geometry. Below we give a description of the limit using the language of the special 
algebra. We start by considering the Kerr-Schild form of the Schwarzschild metric 
written in double null coordinates u = t — x and w = t = x, 

ds^ = dudw - dy"^ - dz^ H kahdx°-dx^ (29) 

r 

where in these coordinates fc" = {{r — x)/r, {r + x)/r,y /r, z/r). 

The Minkowski background enables us to have a well defined concept of boost and 
we may therefore boost the solution by velocity v along the x-axis. We therefore write 



+ l-v_ 
V 1-w \ l + v ^ ' 

and to keep the energy of the "particle" finite we rescale the mass according to the 
special relativistic formula 

m= {l-v^)^p. (31) 

Substituting this into gives a 1-parameter family of metrics depending on the boost 
velocity v. We are interested in the ultrarelativistic limit in which v reaches the speed of 
light (i.e. V — * 1), so we replace t> by 1 — e and regard gab '■= [g{e)ab] as an element of the 
special algebra. It is readily shown that most of the terms in the perturbation ^kakb 
are associated to zero. The only surviving term is ^koko. Using calculations very 
similar to Steinbauer jHI] we find this is associated to 8p\npS{u) and hence gab ~ g{o)ab 
where 

g(o)abdx"'dx^ = 8p\iipS{u)du'^ + dudw — dy"^ — dz^ . (32) 

This metric describes a pp-wave and is fiat everywhere except on the null plane -u = 
which contains the "particle" . This line element was first derived by Aichelburg and Sexl 
P], who started with Schwarzschild written in isotropic coordinates and simultaneously 
boosted the solution and made a w-dependent coordinate transformation to compute 
the limiting metric. 

It should be pointed out that this result is entirely consistent with the calculation of 
the energy-momentum tensor of the Schwarzschild solution given in the previous section. 
The ultrarelativistic limit of the latter in the (u, w, x, y) coordinate system is associated 
to 

5{u)5'^^\y, z)pap' where Pa = (1, 0, 0, 0), (33) 

which is just the energy-momentum tensor of cjab- Indeed this observation was used by 
Balasin and Nachbagauer [Hj to derive the ultrarelativistic limit of the Schwarzschild 
and Kerr geometries (see also [H]). 

It is also possible to calculate the ultrarelativistic limit of the Reissner-Nordstr0m 
solution. However in this case it is also necessary to rescale the charge according to the 
formula 

= (1 - v^)^f, if a constant) (34) 
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in order to obtain a distributional limit at all as v tends to the speed of light. The 
limiting metric was found by Lousto and Sanchez [53 using the methods of Ij to be 



This was confirmed using a calculation in Q by Steinbauer jSUI. The solution obtained 
again represents a pp-wave and is fiat everywhere except on the null plane u = 0. 
The ultrarelativistic limit of the electromagnetic energy-momentum tensor of Reissner- 
Nordstr0m is also found to be compatible with the energy- momentum tensor of ()35|) . 

However the rescaling of the charge has the rather unexpected effect that while the 
electromagnetic field vanishes in the P'-limit the ultrarelativistic energy-momentum 
tensor does not (cf [58 ). Steinbauer in |HI] rephrased this fact using association in Q 
i.e. showed that all the components of the electromagnetic field were associated to zero 
while the 00-component of the energy-momentum tensor was associated with a multiple 
of the delta function. However it must be stressed that regarded as elements of Q the 
electromagnetic field components are non-zero (even though they are associated to zero). 
Calculating the energy-momentum tensor of this field within Q gives the correct result. 
There is nothing unusual within Colombeau theory in having objects in Q which are 
associated to zero having products which are not associated to zero. This is simply a 
reflection of the fact that association does not respect multiplication in general. 

The procedure of Aichelburg and Sexl has also been used to derive the 
ultrarelativistic limit of the Kerr metric by several authors see |57| l59| l2H l88] . In 
fact a number of different sources such as cosmic strings, domain walls and monopoles 
have been boosted to obtain ultrarelativistic spacetimes of impulsive pp-waves which 
in turn have been used to describe (quantum) scattering processes of highly energetic 
particles (see [HZIIHHI for an overview). 

On the other hand Dray and 't Hooft p2l have generalised Penrose's "scissors and 
paste" method [HHl (see also §6) to non-flat backgrounds and used it as an alternative 
way to derive the Aichelburg-Sexl geometry as well as more general gravitational shock 
waves. Using this method Dray and 't Hooft derived the spherical shock wave due to 
a massless particle moving at the speed of light along the horizon of a Schwarzschild 
black hole which was used to study the influence of matter, falling into the black hole 
on its Hawking-radiation. These ideas lie at the heart of 'tHooft's S-matrix approach 
to quantum gravity [86] . 

6. Geodesies for impulsive gravitational wave spacetimes 

In the previous section we showed how the ultrarelativistic limit of the Schwarzschild 
solution lead to an impulsive gravitational wave spacetime. We now consider geodesies 
in such spacetimes with the metric taking the form 




(35) 



ds'^ = f{x^)6{u)du^ + dudw — dAsdx^dx^ 



(36) 
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where A, B = 2, 3 denote the transverse coordinates. These spacetimes have 
been constructed by Penrose using his vivid "scissors and paste" approach (see |7Uj). 
Geodesies for these impulsive gravitational wave spacetimes have been considered by 
Ferrari, Pendenza and Veneziano |23j, Balasin jH] and Steinbauer As one might 

expect, one can regularise the geodesic equations, and show that the geodesies consist 
of broken and refracted straight lines. 

A more rigorous derivation of these results using existence and uniqueness theorems 
within the Colombeau algebra has been given by Kunzinger and Steinbauer in More 
precisely they replaced the delta function in (jHSJ by a generalised function D possessing 
a so called strict delta net {p^} as a representative i.e. they considered the generalised 
line element 



with D = [{pe}] and supp(pe) — ^ {0} and / Pe{x) rfx 1 as e ^ and ||pe||Li uniformly 
bounded in e. They were able to show that the geodesic as well as the geodesic deviation 
equation for the metric (jHTf) may be solved uniquely in Q. Moreover these unique 
generalised solutions possess the physically expected associated distributions. Note that 
diffeomorphism invariance of these results is assured by diffeomorphism invariance of the 
class of strict delta nets. Note further that strictly speaking these calculations have been 
performed using the concept of generalised functions taking values in a manifold, cf §8, 
since geodesies are curves from an interval into spacetime. 

In the literature impulsive pp-waves have frequently been described in different 
coordinates where the metric tensor is actually continuous (cf IZO]). In the special case 
of a plane wave, /(x, y) = — y"^ and m+ denoting the kink function. 



Clearly a transformation relating ()38|) and ()36|) cannot even be continuous, hence in 
addition to involving ill-defined products of distributions it changes the topological 
structure of the manifold. However, the two mathematically distinct spacetimes are 
equivalent from a physical point of view, i.e. the geodesies and the particle motion agree 
on a heuristic level (see 

Using their results on the geodesic equation in Kunzinger and Steinbauer in 
|I3] succeeded in showing that the discontinuous change of coordinates is just the 
associated distributional map of a generalised coordinate transform. More precisely 
modelling the (distributional form of the) impulsive pp-wave metric in a diffeomorphism 
invariant way by the generalised metric (jHTj) the latter may be subject to a generalised 
change of coordinates T. In either coordinates the associated distributional metric 
may be computed giving the distributional (respectively the continuous) form of the 
pp-wave metric. Physically speaking the two forms of the impulsive metric arise as 
the (distributional) limits of a sandwich wave in different coordinate systems. Hence 
impulsive pp-waves are indeed sensibly modelled by the generalised spacetime metric 
fl37|) . However in the different coordinate systems a different distributional picture is 
obtained. 




(37) 




(38) 
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A similar situation arises in the case of impulsive spherical waves which have also 
been introduced in [7D|. In this case the distributional form of the metric arises as an 
impulsive limit of type-N Robinson- Trautman solutions, which however due to the fact 
that the metric is quadratic in the profile function formally involves the square of the 
delta function. An explicit discontinuous coordinate transformation relating this form 
of the metric with the continuous form was given in [721 . ^ study of this situation 
using Colombeau methods rehes on a better understanding of the geodesies of these 
spacetimes; a study of these has been initiated in [7H] . 

7. Conical Singularities and Cosmic Strings 

An important example of a two dimensional singularity is provided by the conical 
spacetime 



where A ^ 1 and is a standard 27r-periodic angular coordinate. This spacetime 
is locally flat for r 7^ and heuristic arguments suggest that it has delta-function like 
curvature on the r = axis. The corresponding energy-momentum tensor then describes 
a string with stress equal to the mass per unit length /i, where /i = 27i{l — A). This is 
precisely the form of the energy-momentum tensor of a cosmic string in the weak field 
thin string limit. Such spacetimes are also important mathematically as they provide 
simple examples of quasi-regular singularities (i.e. singularities for which the components 
of the Riemann tensor measured in a parallely propagated frame tend to a well defined 
limit, see e.g. Vickers lEH! for further details). 

Unfortunately such spacetimes are not gt-regular, so one cannot expect the 
curvature to be well-defined using conventional distribution theory. Furthermore Geroch 
and Traschen showed how it was possible to obtain different values of the mass per unit 
length by taking different regular isat ions. Because of this Clarke, Vickers and Wilson 
[in] chose to investigate the curvature of the cone using the "full" Colombeau algebra 
where one has a canonical embedding of distributions. See §8 below for further details. 

The singular part of the curvature arises from the conical singularity in the 2-cone 



and for simplicity we will describe the calculation of the curvature (density) of this 
metric. Because polar coordinates are not well defined at the origin one must first write 
the metric in a regular coordinate system which includes the origin. To make things 
simple we will choose to work in Cartesian coordinates but as will be shown below the 
final result is independent of the coordinates used. In Cartesian coordinates one has 



ds^ = - dr^ - Ah^dcf" - dz\ 



(39) 



ds^ = dr^ + A^rH<f 



(40) 




(41) 



where 




(42) 
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We now regard the metric as an element of £^a/(R^) by taking the convolution of the 
components Qab with an arbitrary molifier. To find these we need to calculate 



where h{x,y) := e^*"^. 

By expanding in circular harmonics, using the compactness of the support of the 
molifier and using the residue theorem one can obtain a Fourier series for h and hence 
for g(e)ab- One may now estimate the curvature of g(e)ab, and use the Gauss-Bonnet 
theorem to show that 



It is important to stress that the methods described here may be used to calculate 
the curvature of the full four dimensional cone, although one may no longer use the 
Gauss-Bonnet theorem to calculate the curvature and therefore requires more delicate 
estimates (see Wilson jHl] for details). One then obtains the heuristically expected 
energy-momentum tensor. Furthermore it is not hard to modify the results to deal with 
a metric which is not exactly conical, but approaches one quadratically as r ^ 0. A 
second point to note is that the result does not depend upon the coordinates in which the 
calculation is carried out so long as they are smoothly related to Cartesians. This may 
be shown by explicitly transforming to new coordinates X = X{x,y) and Y = Y{x,y) 
and doing the whole calculation in the new coordinates. Again more delicate estimates 
are required but one can show that the resulting curvature density transforms in exactly 
the same way as a delta-function (see Vickers and Wilson |91j| for details). 

The above calculation shows that the curvature of a 4-dimensional cone, when 
calculated in the full Colombeau algebra, has a curvature which is associated to a 
delta-function which gives a mass per unit length equal to the deficit angle as one 
would expect from the heuristic calculation. Note that this result does not say that the 
curvature is equal to a delta function, but simply that it is associated to a delta function. 
Thus if one works at the level of association this result shows that regularisations based 
upon smoothing convolutions give an unambiguous answer for the curvature of the cone. 
Indeed a number of authors have applied different regularisations and obtained the same 
result. Balasin and Nachbagauer |lj used a regularisation based on a family of smooth 
hyperboloidal surfaces converging to a cone. A number of authors including Marder 
|?)T] and Geroch and Traschen [22] have looked at "rounding off" the point of a cone 
with a spherical cap, and Louko and Sorkin jSH] have looked at a regularisation based 
on a different coordinate system. Unfortunately things are not quite as straightforward 
as one would hope. It is also possible to choose regularisations which do not yield a 
mass per unit length equal to the deficit angle. An important example of this was given 
by Geroch and Traschen [22] • They first formed a regularisation sequence g{e)ab which 
gave the standard answer for the mass per unit length. They then introduced a second 
regularisation sequence g(e)ab which was related to the first sequence by a conformal 




(43) 



[R^/^] ^47r(l - A)5(2)_ 



(44) 
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factor 



g{e)ab = g{e)ab, 



(45) 



where 



n = exp(A/(r/e)) 



(46) 



and / is a smooth function whose support is [1/2, 1]. They then found that the mass 
per unit length of the hmiting spacetime was dependent on both A and /: 



From the point of view of Colombeau theory this "bad" behaviour of the above 
regularisation is not unexpected since although fi^ ^ 1, it is not equal to the unity 
function in ^(M^) and hence does not represent the conical spacetime in the full 
Colombeau algebra. Wilson jHlj looked at general regularisations of conical spacetimes 
and gave a condition that ensured that a given regularisation sequence gave the standard 
result for the mass per unit length. He showed that provided that the difference between 
the connection one-forms of the regularised spacetime and the original spacetime 
diverged no faster than 1/r then the answer would be the standard result. Thus 
although it is possible to obtain regularisations which give a different mass per unit 
length, the geometry of these bad regularisations diverges strongly from that of a cone 
as one approaches the axis. 

A significant generalisation of the curvature calculations for a conical spacetime 
was obtained by Wilson [OH] who extended the results to a four dimensional time 
dependent cosmic string. He considered cylindrically symmetric pure radiation solutions 
of Einstein's equations. These metrics may be written in the form 



and naively one would expect a delta-function contribution to the curvature due to 
the angular deficit of 27r(l — e"'^^*-'). However, because of the time dependence of the 
angular deficit, the singularity cannot be quasi-regular but must be stronger. In fact 
it is an example of an "intermediate singularity" and the components of the Riemann 
tensor have a limit in a special choice of frame. Wilson calculated the energy-momentum 
tensor of this spacetime by first writing the metric in null Cartesian coordinates and 
then proceeded to obtain a smooth metric by convolution, estimate the curvature of the 
smooth metric and show that it is associated to a distribution. This result shows that 
even in the radiating case the mass per unit length is given by the expected formula 



Furthermore by applying the methods of Ashtekar et al. [S] one may show that this 
result agrees with the asymptotically measured mass per unit length. However unlike 
the case of the static string it is unclear whether the calculation is coordinate invariant. 
The (m, X, y, z) coordinates used are natural for investigating radiating systems but have 
the disadvantage that Minkowski space appears to be singular on the axis when written 




1/2 



(47) 




(48) 




(49) 
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in these coordinates. However the correct mass per unit length for a static cosmic string 
is obtained, and this vanishes for Minkowski space. 

8. Nonlinear distributional geometry 

An underlying principle of general relativity is that the measurement of physical 
quantities should be independent of the coordinate system used. Mathematically this is 
reflected in the fact that the theory is formulated in terms of tensor fields on manifolds. 
The minimum differentiability of the coordinate transformations and the dependence on 
the differential structure of the manifold is quite subtle, but at the very least the theory 
should be invariant under smooth diffeomorphisms. 

The definition of the special algebra introduced in §2 may be generalised in a 
more or less straight forward manner to yield a special Colombeau algebra Q{M) on a 
differentiable manifold M [21]. This setup was extended by Kunzinger and Steinbauer 
in to a theory of generalised sections of vector bundles. Furthermore in the 
study of (pseudo-)Riemannian geometry in the generalised setting was initiated and in 

generalised connections and curvature in general principal and vector bundles were 
studied. However both geodesies and diffeomorphisms involve considering functions 
with values in a manifold and if one wishes to consider these objects in the generalised 
setting one is forced to consider manifold valued generalised functions (a concept which 
it is not possible to deal with using distributions). The study of such functions was 
first looked at in |1T], where the space ^[X, F] of generalised functions defined on the 
manifold X and taking values in the manifold Y was defined. This work was extended 
to a functorial theory in |17] where several global characterisations of the notions of 
moderateness and negligibility for generalised functions from X to y are given. These 
characterisations provide the key to proving that composition of generalised functions 
between manifolds can be carried out unrestrictedly. Using these ideas it is possible to 
consider generalised geodesies and flows as well as singular ODEs on manifolds (see |lHj 
for details). An overview over this "nonlinear distributional geometry" can be found in 
[2H1 Ch. 3.2]. 

However in all these constructions the embedding of distributions and functions 
of finite differentiability is non-canonical; in addition to the dependence on the choice 
of a molifier (see §2), any embedding into Q{M) will not be diffeomorphism invariant. 
So how can this setting be of use in a diffeomorphism invariant theory like general 
relativity? 

Firstly in some applications there will be a natural physical parameter, such as a 
coupling constant, that may be used to directly provide a representation of the singular 
objects involved in Q. Hence there will be no need for using any embedding at all. 

Another possibility is to model a singular spacetime metric by a whole class of 
generalised metrics that is by definition diffeomorphism invariant. This has been done 
in case of the description of impulsive pp waves in Q (cf §6). 

Finally one may regard the Colombeau algebra as an intermediate calculational 
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tool for obtaining distributional answers. One picks some coordinate system, embeds 
the components of the metric into the Colombeau algebra (in the given coordinates) 
and calculates the components of the curvature and energy-momentum tensor etc. One 
then tries to show that these objects are associated to distributions. Finally one repeats 
the entire calculation in some other (smoothly related) coordinate system and tries to 
show that the answer transforms in the expected way for a distribution. 

Although such calculations can be useful in some situations it would be preferable 
to make the embedding into the algebra coordinate invariant. To explain how this can 
be done we first introduce the full Colombeau algebra Q^{W^) which, unlike the 
special version, has a canonical embedding of distributions. 

This is achieved by substituting the index set (0, 1] by a suitable class of molifiers. 
More precisely we introduce the following grading on the space of all molifiers 

Ao := {$ G ViJT) : j <l>(x) dx = 1} 

Ag := {$ G A : J $(a;)x" dx = , 1 < \a\ < q} (g G N) 

and take the basic space to be 

■= {/ : A X M" ^ M"| / smooth in x}. 

Now the respective spaces of moderate and negligible functions may be defined as 
follows (again ^^(a;) := (l/e"')$(a;/e)). 

8.1 Definition. 

(i) ( Moderate functions ) 

£^(M") := {f eS^: 'iK CC W "ia G 3p G Nq G ^ : 

sup |D"/(<l>„x)| = 0{e-P) as 6^0} 
xeK 

(a) ( Negligible functions ) 

A/'"(M") := {/ G S{n) : \/K CC M" Va G N[} Vp G No 3g V<l> G A : 
sup |D°i?(<l>e,x)| = 0(eP) ase^O} 

(Hi) (Full algebra) 

Distributions are now simply embedded into by convolution with the molifiers i.e. 

i{T) = [T * $] (50) 

and one obtains a differential algebra of generalised functions on M" (or open subsets 
thereof) just as in the special version with the additional benefit of a canonical 
embedding of the space of distributions. 

Unfortunately this construction cannot be generalised to the manifold setting in 
a simple way as the definition of the spaces Aq is not invariant and moreover the 
embedding is not diffeomorphism invariant since convolution again depends on the linear 
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structure of M". However an invariant embedding can be achieved by demanding that the 
mohfiers $ transform in an appropriate way. Colombeau and Meril in made the first 
decisive steps towards a diffeomorphism invariant full Colombeau algebra by weakening 
the moment conditions to only hold asymptotically (which makes them invariant) and 
enlarging Aq to a space of bounded paths e i— >• G P(]R"). A flaw in this construction 
was found and removed by Jelinek |3H] who developed an improved version of the theory 
which involved some subtle changes of definitions and established a number of important 
technical results. These ideas were then fully developed and given a firm mathematical 
basis in two papers dealing with the foundations of nonlinear generalised functions fI7\ 
where the first diffeomorphism invariant full Colombeau algebra on (open sets of) M"^ 
was constructed. 

For applications in general relativity it is moreover desirable to have a geometric 
and global version of the theory rather than simply giving transformation rules for the 
local theory. Such a construction was given in [211 (for an overview see ISO])- Here we 
only mention that the key idea is to replace the scaled and unbounded paths (l/e)$'' 
which are employed in the definition of moderateness and negligibility in the local theory 
by smoothing kernels $ which are C°° maps from (0, 1] x M to compactly supported 
n-forms on M. In this way one obtains a geometrically constructed full Colombeau 
algebra on a differential manifold M where the canonical embedding of distributions 
commutes with Lie derivatives. 

However this theory still lacks a canonical embedding of distributional tensors. 
Although the work of jSHI described above provides one with an invariant embedding 
of scalars, one cannot simply apply this to the components of a tensor and obtain 
an invariant embedding of the tensor. Embedding the components of a tensor and 
then transforming would in general give a different answer from transforming and then 
embedding since multiplication by a smooth function does not in general commute with 
the embedding. The problem really stems from the way the convolution integrates the 
components of the tensor at different points of the manifold. The solution to this is 
to introduce some additional structure which enables one to first transport the tensor 
fields to the same point p in M so that one can then do the integration in a meaningful 
way. Such a transport is naturally provided by specifying a background connection, and 
in keeping with the spirit of the full algebra this is made an argument of the generalised 
tensor field. Thus a generalised tensor field depends upon a background connection 7, 
a smoothing kernel $ and the point p. This is currently work in progress but see [OHj 
and 35] for a more detailed description of this approach. 

This approach now provides a "nonlinear distributional geometry" plus a canonical 
and invariant embedding of distributional objects. For example one can canonically 
embed the conical metric of the previous section to obtain a generalised metric. One 
then finds that this metric has a curvature tensor which is associated to the Dirac 
distribution and this shows that in any coordinate system a conical spacetime has a 
generalised Einstein tensor which satisfies 

Gab ~ fab (51) 
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where Tab is the embedding into the Colombeau algebra of the energy-momentum tensor 
of a (thin) cosmic string. 

9. Weak singularities and Generalised hyperbolicity 

As discussed earlier the definition of a singularity in general relativity is different from 
other field theories where one has a background metric. In general relativity one 
detects the presence of singularities by showing that the spacetime is incomplete in 
some sense. In the standard approach to singularities (see e.g. Hawking and Ellis, 
|H2]), a singularity is regarded as an obstruction to extending a geodesic. However this 
definition does not correspond very closely to ones physical intuition of a singularity. 
This led to a consideration of whether physical objects would be subjected to unbounded 
deformations as one approached the singularity and was formulated mathematically in 
terms of strong curvature conditions. Unfortunately it is hard to model the behaviour 
of real physical objects in a strong gravitational field and because of this Clarke [T^ 
suggested that one should consider instead the behaviour of physical fields (for which 
one has a precise mathematical description) near the singularity. According to the 
philosophy of "generalised hyperbolicity" one should regard singularities as obstructions 
to the Cauchy development of these fields rather than an obstruction to the extension 
of geodesies. However even a mild singularity is an obstruction if one uses the standard 
theory of distributions. By considering solutions to the wave equation in the Colombeau 
algebra one finds a certain class of weak singularities which do not prevent the evolution 
of test fields (Vickers and Wilson, 93 ). 

In this section we will look at solving the wave equation on a spacetime with a 
locally bounded singular metric (such as the conical spacetime considered in §7). The 
standard Cauchy problem takes the form 

□M(t,x") = (52) 
m(0,x") =t;(x") 
9iM(0,x") = w;(x") 

with initial data {v,w) lying in the Sobolev spaces H^{S) x H^(S) prescribed on the 
initial surface S, given by t = 0. Because of the form of the metric one would expect the 
solution (if it exists) to be defined as a distribution. This however will cause difficulty 
in interpreting 

Du = {-g)-y^da{{-gY'^g'^'dbu) (53) 

as a distribution in the framework of classical distribution theory because the above 
equation has (non-constant) singular coefficients and involves ill-defined products. Again 
we overcome these difficulties by using the nonlinear generalised function theory of 
Colombeau. For an overview of the treatment of PDEs with singular coefficients, data 
and solutions in this setting see [B^. 

We first canonically embed the metric Qab into the full Colombeau algebra Q^{M) 
by using a convolution integral (fTHjl as in §7 to obtain a representative {cj(e)ab) £ ^ai{M). 
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Since the initial data {v, w) does not have to be smooth, we must also embed it into 
the algebra as (V, W) represented again by convolution integrals denoted and 
respectively. 

The generalised function wave operator acting on a generalised function U 
represented by may then be written as 

CTw, = {-g,)-'/'da{{-g.Y/'g:%u,). (54) 

We would like to then be able to solve the Cauchy problem in the space Q^{M) 

□t/(t,x") = (55) 
[/(0,x") =F(x") 
dtU{0,x'') = Vr(x") 

and obtain a solution U G Q^{M) which is associated to a distribution. In practice one 
works with the equivalent problem in ^^(M) 

CTn,(t,x") = /,(t,x") (56) 
M,(0,x") =f,(a;") 
9iM,(0,a;°) =w,(x") 

where (/J is negligible. 

In inn] it is shown how to estimate solutions of (jKHj) and its derivatives in terms 
of powers of e, given the moderate and null bounds of /e, and w^, using a method 
of energy estimates following Hawking and Ellis [22j| and Clarke The additional 

complication in this situation is that one needs more explicit bounds because one needs 
to know the precise way in which the constants depend upon e. This is accomplished by 
working with function spaces defined using higher order energy estimates (related to the 
super-energy tensors of Senovilla [TH]) for which one has bounds in terms of the covariant 
derivatives of the curvature. Using this method one can show that the Cauchy problem 
()55|) has a unique solution U G Q^{M) which is independent of the representation chosen 
for gab, V or W. In the case of a conical spacetime one can go further and show that 
this solution is actually associated to a distributional solution. Thus a conical spacetime 
satisfies the condition of "generalised hyperbolicity" as claimed. It seems likely that 
these techniques may then be used to show that a much wider class of spacetimes with 
weak singularities satisfy the generalised hyperbolicity condition. See ([H3j) for some 
recent results in this direction. 

10. Conclusion 

In this review we have looked at the extent to which it is possible to use conventional 
distribution theory to look at solutions of Einstein's equations. Although there is an 
important class of new solutions that can be obtained by going beyond the confines 
of C^" metrics the largest class that one can work with that is "stable" is given by 
the gt-regular metrics. Such metrics can be used to describe solutions with singular 
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support of the curvature on a hypersurface but are unable to deal with singularities 
of higher codimension. To deal with these it is necessary to go beyond distributions 
and work with a theory of nonlinear generalised functions. We have shown that such 
an appropriate description of nonlinear generalised functions is given by the theory of 
Colombeau algebras. 

The special theory provides a straightforward computational tool for calculating 
the distributional curvature of a number of singular metrics and throwing some light 
on the physical nature of the singularity. Furthermore it has also been possible to 
define generalised functions taking values in a manifold and this allows one to talk 
about generalised geodesies and generalised symmetries (see e.g. 0) of a spacetime. 
Unfortunately due to lack of space we have had to omit from this review both this 
latter topic and the topic of impulsive pp-waves and ultrarelativistic black holes with 
non- vanishing cosmological constant (see e.g. [71J). 

However the special algebra does not provide one with a canonical embedding of 
distributions so there is always a question about the extent to which the answer depends 
upon the particular embedding that is used. The full Colombeau algebra rectifies this 
problem and a global formulation that is independent of the coordinate system has 
been given. Although the details of the tensorial theory remain to be fully worked out 
this work provides the basis for a coordinate and embedding independent theory of 
generalised (pseudo-)Riemannian geometry which can be used to analyse a wide class 
of singular spacetimes. In particular it is possible to give a distributional interpretation 
to a many physically reasonable singularities. The remaining singularities can therefore 
be regarded as true gravitational singularities. An outstanding project is to consider 
the singularity theorems in this generalised setting and show that they predict true 
gravitational singularities rather than simply distributional singularities. 
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